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I. INTRODUCTION 


Unlike ground-based gravitational wave (GW) detectors f4]-[6], space-borne GW detectors [DIMS] usually have 
varying interspacecraft distance and the clock based in each spacecraft (S/C) [as the vertex of the detector] are 
freely running individually. To achieve the extraordinary target sensitivities of the detectors, precise determination of 
interspacecraft distances and accurate synchronization of the separate spacecraft clocks in different S/C are usually 
necessary. This happens in the first stage of data processing on ground. 

We have tried to solve this for (e)LISA pfl [7H9] in the preprocessing stage in our previous paper [3], where a 
hybrid-extended Kalman filter was designed to preprocess raw interspacecraft measurements of (e)LISA. Although 
this ranging determination and clock synchronization problem shares some similarity with the global positioning 
system (GPS) track problem [16, .17:, the differences [3] between them make it far more complicated than the GPS 
case, see [3] for detailed discussions. The 24-dimensional state vector designed in [3 j cannot be determined well due 
to the insufficient number of measurements. In this paper, we will explore alternative dynamic models for Kalman 
filters to reduce the ambiguity in the state vector, while using the same measurements. We investigate in detail how 
different filtering dynamic models affect the accuracy of the estimates. These various dynamic models also provide 
better means to deal with other effects, such as jittering recording time, or the Sagnac differential delay. 

In Sec. |TTJ we describe the inter-spacecraft measurements that are going to be used to improve the ranging accuracy 
and to synchronize the clocks in dif fere nt S/C. We also review the full 24-dimensional physical state vector used for 
Kalman filters in paper [3]. In Sec. Ill and |IV| we try to model the system with 23-dimensional and 22-dimensional 
state vectors by excluding the common absolute time errors and the frequency errors of the clocks. In Sec. [Vj we 
describe a toy system model with only the clock variables. Although this toy model is too simplified for realistic 
simulation, it allows us to see the analytic forms of all the components of the measurement matrix H and the dynamic 
matrix F. It may provide insights to more complicated problem in reality. In Sec. |VI| and |VII[ we explore two 
effective models with phenomenological parameters to further reduce the ambiguity in the state vector. The effective 
models try to use variables that are directly measured, thus reducing the nonlinearities in the measurement equations 
and increasing the numerical accuracies of the Kalman filter algorithms. The number of variables in the state vector 
of these effective models is also significantly smaller, and their matrices are better conditioned. It turns out that the 
effective models are simpler and more efficient. Finally, we summarize the results in Sec. VIII 


II. INTERSPACECRAFT MEASUREMENTS AND THE FULL STATE VECTOR 


The interspacecraft measurements o Dam] of (e)LISA contain information about the absolute distance between 
the S/C, and their clock errors. This information is expected to be extracted from the measurements in the pre¬ 
processing stage [3] of (e)LISA data analysis. At S/C j, the local carrier laser with frequency J9 arrier beats with the 
incoming laser from S/C i. The beatnotes Dij contain information of the Doppler shifts, hence we call them the 
Doppler measurements, which are given as follows 


Dij — 


^carrier 


^carrier 


+ /V1 


1 - 


1 - 


(Vj -Vi) 


Sfj 

Jnom 




(i) 


where is the velocity of S/C i, c is the speed of light, /9 arrier denotes the carrier frequency of the laser from S/C i, 
fiij denotes the unit vector along the axis from S/C i to j, /f om and Sfi are the nominal frequency and the frequency 
error of the clock, and nR is the measurement noise. Since this measurement contains the expected gravitational wave 
signals /^ w , it is also referred to as the science measurement. 

Pseudo-random codes are modulated onto the carrier lasers to measure the distance between the S/C [20] . By cor¬ 
relating the remote laser with an electronic copy of the same pseudo-random code in the phasemeter, the propagation 
time of the laser is determined, similar to the pseudorange in the global navigation satellite system [T6l H7] . We call 
these the ranging measurements and denote them as follows 

Rij = Lij + (STj — STi) • c + nfj, (2) 

where denotes the arm lengths, and STi denotes the time error of the clock in S/C i, which is directly related to 
the frequency error [3|. 

The clock frequencies are up-converted to the GHz range and also modulated onto the carrier lasers. The sidebands 
thus created are called the clock sidebands. They contain the information of the frequency errors of the clocks. They 
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measure the same frequency change as the carrier-to-carrier beatnotes, and, in addition, they measure the amplified 
differential clock errors. Here, we use the simplified formulae of the clock measurements obtained after subtracting 
the common-mode component Dij and dividing by the up-conversion factor [3] 

Cij = Sfj — Sfi + nfj. (3) 

More detailed descriptions of the interspacecraft measurements can be found in [3] and the references therein. 

By writing the subscripts explicitly, we have 18 measurements in total 

V = (-^ 31 ? Csi , i?2ij ^21, C21, ^12, D 12, C12, 

^32, ^32, C32, i?23> ^23, ^23, ^13? ^13? Cis ) T . 

However, to fully describe the system in this aspect by physical variables, we need a 24-dimensional column state 
vector |3| 


x = (x 1 ,x 2 ,X 3 ,v 1 ,v 2 ,v 3 , 6 T 1 ,ST 2 , 6 T 3 ,Sf 1 , 6 f 2 ,Sf 3 ) T , 

where x% = (xi,yi,Zi) T are the positions of the S/C in the solar system, and Vi = {v x i,v y i,v z i) T their velocities. It 
is clear at this point that we have more unknown variables than measurements in the Kalman filter model. This is 
obvious from the fact that the relative measurements between the S/C do not provide information on the absolute 
positions and velocities of the S/C. 


III. A KALMAN FILTER MODEL WITH A 23 DIMENSIONAL STATE VECTOR 


In particular, the clock errors 5Ti appear in the measurement equations only in the form of differential time errors 
STj — STi. Therefore, the commom mode of the clock errors (8T\ + 8T2 + 5 Ts )/3 cannot be determined. To eliminate 
this degeneracy in the Kalman filter model, we replace the three clock errors <STi, ST 2 and 5T^ by two differential 
clock errors 5 T\ — 8T2 and 8T2 — ST3. In consequence, the state vector is now a 23-dimensional column vector 


x = (fi,f 2 ,^3,^i,^2,^ 3 ,^Ti - 5T 2 ,5T 2 - ST 3 , 

Sfi,Sf 2 ,6f 3 ) T . 

The third differential clock-error variable 8 T% — STs can be expressed as ( 8 T\ — 8 T 2 ) + ( 8 T 2 — 8 T 3 ). The dynamic 
equation for the clock errors can simply be modified as 

— (ST- — 6 T-) = _ (4) 

l ) j nom j nom ' v / 

The dynamic matrix and the observation matrix Hk are then modified accordingly. 

We carried out simulations to compare the performance of this Kalman filter model with the performance of the 
model designed in paper [3] with a 24-dimensional state vector. The two filters have been run over the same simulated 
measurement data. Fig.[l]shows a comparison of the arm length determination between the two Kalman filter models. 
It can be seen that the performance of the two Kalman filter models is comparable in determining the arm lengths. 
Both models can successfully decouple the arm lengths from the initial clock biases and reduce the noise in the arm- 
length roughly by one order of magnitude. Notice that the performance of a Kalman filter depends on the specific 
noise realization. Therefore, a small difference between the two models in the estimation error is insignificant. Fig. [2] 
shows histograms of estimation errors in differential clock errors. Both models are able to reduce the differential clock 
errors roughly by one order of magnitude. The Kalman filter model with a 23-dimensional state vector performs 
slightly better. Fig. [3] shows histograms of estimation errors in differential frequency errors. Both Kalman filter 
models have greatly reduced the noise in the raw measurements. The performances turn out to be comparable. 

All in all, the Kalman filter model with a 23-dimensional state vector designed in this section performs slightly 
better than the model with a 24-dimensional state vector. The reason is that the Kalman filter model with a smaller 
state vector has reduced the ambiguity in the system model, while retaining the full information on the measurement 
mechanism. 
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arm length precision [m] arm length precision [sec] 


(a) (b) 

FIG. 1: A comparison of the arm length determination. Fig. (a) shows histograms of errors in raw arm-length 
measurements and Kalman filter estimates with a 24-dimensional state vector. Fig. (b) shows histograms of errors in 
raw arm-length measurements and Kalman filter estimates with a 23-dimensional state vector. Notice that the 
initial clock biases are not included in the raw measurement errors for better vision. 




relative clock jitters [sec] xIO’ 9 relative clock jitters [sec] 


(a) (b) 

FIG. 2: Histograms of estimation errors in differential clock errors. Fig. (a) shows histograms of errors in the raw 
data and Kalman filter estimates with a 24-dimensional state vector. Fig. (b) shows histograms of errors in the raw 
data and Kalman filter estimates with a 23-dimensional state vector. 


IV. A KALMAN FILTER MODEL WITH A 22 DIMENSIONAL STATE VECTOR 

In this section, we try to further simplify the state vector. From the simulation, we know that the absolute frequency 
error Sfi cannot be determined accurately, since the common mode of the three frequency-error variables is poorly 
determined by the measurements. (Here, the ‘absolute’ frequency error is in contrast with the differential frequency 
error Sfi — Sfj .) The only information about the absolute frequency errors comes from the term 1 — Sfj/ff om of Eq. 
|TJ Usually, Sfj/ff om is several orders of magnitude smaller than 1, thus this term provides very limited information. 

We can actually approximate this term by 1 without losing much information. Then, the Doppler measurement 
equation becomes 

^i) + fS w 

+»g. (5) 

Now, in all the three interspacecraft measurement equations [2j [5] and [3] the frequency errors only appear in the form 
of differential frequency errors Sfj — Sfi. Therefore, we can replace the three frequency error variables Sfi by two 
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frequency [Hz] frequency [Hz] 

(a) (b) 

FIG. 3: Histograms of estimation errors in differential frequency errors of the clocks. Fig. (a) shows histograms of 
errors in the raw data and Kalman filter estimates with a 24-dimensional state vector. Fig. (b) shows histograms of 
errors in the raw data and Kalman filter estimates with a 23-dimensional state vector. 



histogram of frequency difference 



(a) 


(b) 


(c) 


FIG. 4: Simulation results of the Kalman filter model with a 22-dimensional state vector. Fig. (a) histograms of 
errors in raw arm-length measurements and Kalman filter estimates. Fig. (b) shows histograms of the raw and 
residual differential clock errors. Fig. (c) shows histograms of the raw and residual differential frequency errors. 


differential frequency error variables Sfi — 6/2 and Sf 2 — Sfo. The state vector is correspondingly reduced to 22 
dimensions: 


x = (x 1 ,x 2 ,x 3 ,v 1 ,V2,V3,5T 1 - ST 2 , 6 T 2 - ST 3 , 
Sfi-Sf 2 ,Sf 2 -Sf 3 ) T . 


The dynamic model for the clock errors and the frequency errors is approximated as follows 


d 

~5Tj - 6Ti 


'2 {sfj - +fr m ) 

d t 

Sfj - 5fi_ 


0 


The main advantages of this Kalman filter model are (i) the further reduction of the near degeneracy in the system 
model, (ii) the reduction of the nonlinearity in the measurement equations. On the other hand, the use of approximate 
measurement equations could also be disadvantageous at the same time. 

We implemented this Kalman filter model to process the same simulated measurement data used in the last section. 
We have set the nominal frequencies of the clocks in different S/C to be identical in the simulation, i.e. f^ om = /^ om . 
If the nominal frequencies are not identical, the performance of this Kalman filter model would be degraded depending 
on how different the nominal frequencies are. The simulation results are summarized in Fig. [4j where the estimation 
errors in the arm lengths, clock errors and their frequency errors are shown. Comparing these results with Fig. M 
and [3j we find that the estimation errors of this model in the arm lengths are slight larger than those of the other 







































6 


two Kalman filter models. However, this Kalman filter model performs slightly better than the other two models in 
determining the clock errors. Overall, the performances are similar. 

V. A TOY KALMAN FILTER MODEL WITH ONLY CLOCK VARIABLES 


In this section, we show a toy Kalman filter model only for the clocks. It is simple and can be implemented efficiently. 
However, this model is accurate only when the interspacecraft laser links are instantaneous, hence not straightforward 
to be generalized to the complete treatment with relativistic interspacecraft delays. At any rate, in this model we 
need not to deal with large matrices with hundreds of nonlinear components but only small constant matrices. The 
benefits of studying such a toy model are the following: (i) the matrices in this model are of much smaller size and can 
be shown explicitly with clear dependence on the variables, hence it helps the audience to understand how Kalman 
filters work for (e)LISA; (ii) with the clock variables separated from the arm-length and the Doppler measurements, 
it can be used to design and tune sophisticated clock models more efficiently. 

We first define a 4-dimensional state vector 


x = 


~STi - <5T 2 - 


~Xi 

ST 2 - 5T 3 


%2 

Sfi ~ 


%3 

S' 

1 

On 

co 


_X 4_ 


(7) 


which consists of only the differential clock errors and the differential frequency errors. The dynamic equations can 
be simply modeled as 



~xi 


fx 3 /f nom l 

d 

x 2 


X4 // n ° m 

d t 



0 


_X 4 _ 


0 


noise, 


where we have assumed the nominal frequencies of the three clocks are identical, 
dynamic equations are linear, the dynamic matrix is simply 


i.e. ff 


(8) 


Since these 


F = = 


5/ 

dx 


■0 0 l// n( 
0 0 0 
0 0 0 
0 0 0 


0 

V/ n( 

o 

0 


(9) 


Next, we set up the measurement equations, 
obtained in the following way 


For instantaneous laser links, the differential clock errors can be 


STj STi — :^ c (Rij Rji 


( 10 ) 


where the ranging measurements Rij are given in Eq.j2l The differential frequency errors of the clocks are indirectly 
inferred by the modified clock measurements, see Eq. [3f As mentioned before, the Doppler measurements (see Eq. [l]) 
contain little information on the frequency errors, hence for this Kalman filter model with only clock variables they 
are almost irrelevant. Therefore, the measurement equations are simplified to the following linear form 


V = 


which leads to a constant measurement matrix 


-< 5 Ti - ST 2 - 


Xi 

5T 2 - 5T 3 


X2 

6 T 3 - 


1 

1 

to 

5fi ~ Sf 2 

+ noise = 

X3 

Sf 2 - Sf 3 


£4 

Uf 3 ~ Sf\ J 


-x 3 - X 4 _ 


h =T 

OX 


-1 -1 0 
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noise, 


(ii) 


( 12 ) 
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We implement the above simplified Kalman filter model to process the same simulated measurement data used in 
previous sections. The simulation results are shown in Fig. [5] It can be seen that in the considered case the simplified 
Kalman filter has significantly reduced the noise in the differential clock errors and greatly reduced the noise in the 
differential frequency errors. Comparing these results with Fig. @0 and [4] the performance of this simplified Kalman 
filter is a few times worse than other Kalman filter models in determining the relative clock jitter and about ten times 
worse in determining the differential frequency errors of the clocks. 


VI. A PERIODIC SYSTEM MODEL 

If each S/C of LISA would follow a Kepler orbit, the variations of the arm lengths of LISA were strictly periodic. 
Even for relativistic orbits with gravity contributions from all the planets in the solar system, the variations of the arm 
lengths are also mostly periodic. Using the Kepler orbital setup in the paper [21], the annual evolutions of the LISA 
arm lengths are shown in Fig. [6] (a). The arm-length variations resemble sinusoidal functions within 6 months. In the 
subsequent half year, the arm length variation remains a roughly sinusoidal shape, but with a different amplitude. 
By optimizing the inclination angle between the ecliptic plane and the LISA constellation plane [22], the arm-length 
variations can be further reduced, and we obtain the arm-length variations shown in Fig. [6] (b). These arm-length 
variations are similar to a sinusoidal function with a period of one year. In either cases, we can phenomenologically 
model the arm-length evolutions as follows 


L(t) = L + Asin(co’t), (13) 

where L is the average arm length, and uj is angular frequency with either a half-year period or a one-year period. If 
we define the arm length change as below 


A L(t) = L(t) — L = Asin(u;£), 


the dynamic equation of the arm length change is the same as that of a simple harmonic oscillation 


d 2 A L 
d t 2 


T cj 2 AT — 0. 


(14) 


(15) 




relative clock jitters [sec] xIO ' 9 frequency [Hz] 

(a) (b) 

FIG. 5: Simulation results of the simplified Kalman filter model with only clock variables. Fig. (a) shows 
histograms of errors in the differential clock errors. Fig. (b) shows histograms of errors in the differential frequency 

errors of the clocks. 


A Kalman filter based on this simple periodic dynamic model can thus be designed. We define a 11-dimensional 
state vector as follows 


X — (AL 2 i, U 2 l, AT 32 , U 32 , ATi 3 , U13, 

- ST 2 ,5T 2 - 5T 3 , Sh,6f 2 ,Sf 3 ) T . 






















(a) 


(b) 


FIG. 6: The annual arm-length variations of the LISA constellation for Kepler orbits. Fig. (a) has fixed the 
inclination angle between the ecliptic plane and the constellation plane of LISA as 60°. Fig. (b) has varied and 
optimized this inclination angle to minimize the annual arm-length variations. 

The dynamic equations for the arm-length variations and the relative tangential velocities can be obtained from Eq. [15] 
by rewriting it as first order differential equations 


d 
d t 

The dynamic equations for the clock variables are given as follows 
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noise. 


(17) 

In order to make the system model fit better into the Kalman filter frame, we subtract the average arm length from 
the ranging measurements, thus obtaining the following observation vector 

y' = (R31 — A, D 3 i, C31, R21 — £, D 2 i, C21, R12 — A, 

D 12, C12, R32 — T, T> 32 , C32, R23 — L, 

R*23i C 2 3, R\3 — Z, -D 13 , Gi 3 ) T , 


where the interspacecraft measurements are 

Rij = Lij + (STj — STi)c # nfj, 

^carrier _ ycarrier ^ ^ _|_ j 


Dij — 


GW 

ij 


1 - 


Sfj 

j nom 




D 
ij ’ 


Cij = Sfj - Sfi + nF. 


The Measurement matrix H' k = 


is 18-by-ll. We explicitly give its first three rows below 
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where we have omitted the step index k in the components of the matrix. Notice that the dynamic matrix for 
this model is a constant matrix, and most components of the measurement matrix H' k are constant. The nonlinearity 
of the Kalman filter model has been significantly reduced. Also, the number of the variables in the state vector is 
11, which is less than the number of measured quantities 18. Except for the absolute frequency errors of the clocks, 
all other variables in the state vector are directly constrained by the measurements, hence in principle they can be 
determined by the Kalman filter. 

We design a hybrid-extended Kalman filter based on the model described above to process simulated LISA measure¬ 
ment data. The estimation errors in the differential clock errors, arm lengths and the differential frequency errors are 
plotted in Fig. [7] With this phenomenological periodic model, the designed Kalman filter has successfully estimated 
these three kinds of interspacecraft quantities. We also find that the designed Kalman filter model is not sensitive to 
the actual orbits of the LISA constellation. With either angular velocities u) = 2tt rad/year or u) = 4tt rad/year and 
with either LISA orbits shown in Fig. [6] (a) or (b), the performance of the designed Kalman filter turns out to be 
similar in terms of estimation errors. 


VII. AN EFFECTIVE SYSTEM MODEL 

In this section, we wish to design another effective system model, whose errors can be accessed analytically. Com¬ 
pared to the periodic system model, this form of system model is expected to be directly applicable to relativistic LISA 
orbits under consideration of all the planets in the solar system. The arm lengths of LISA are smooth, slowly varying 
and (nearly) periodic functions of time, no matter whether they are calculated in the non-relativistic or relativistic 
framework. The smooth arm-length functions of time can be decomposed into harmonics 

oo 

L(t) = L + An sin (nut + (f> n ), (21) 

n= 1 

where L is the average arm length, and A\ is the lowest order arm-length variation, whose value is about Z/100 
according to the orbit design of LISA. 

The arm-length function can be expanded in polynomials around any time to 

oo oo 

L(t) = L + A n sin(nci;to + (j) n ) + A n cos(nu;to + (j) n )(nujAt) 

n=1 n=1 

^ oo 

-- y A n sin (nuit 0 + 4> n )(nuAt) 2 + 0[Ai 3 ], (22) 

n= 1 

where At = t — to- For elliptical Kepler orbits, A n decays exponentially with n. For relativistic orbits, A n also 
decays much faster than linearly, which leads to (n + \)A n +% nA n . For a short time At = 1000s, we estimate the 
contribution of each order as follows 


A\ujAt rsj 10 4 m, 

(23) 

^Ai(ujAt ) 2 ~ lm, 

(24) 

^Ai(a;At) 3 ~ 10 -4 m. 

(25) 


As a consequence, if we want to design a Kalman filter to process 1000 s LISA measurement data, the following 
polynomial model characterize the LISA arm lengths to 0.1mm accuracy. 

oo 

L(t) = L + A n sin(ncjto + </> n ) + vAt 

n= 1 

+ + 0[At 3 ], (26) 

where v and a are phenomenological variables. 

According to the above phenomenological model, we define a 14-dimensional state vector in the Kalman filter as 
follow 


X — (1/21, ^21, «21, A 32 , ^32, ^32, ^13, rq 3 , a 13 , 

ST! - ST 2 , 5T 2 - 5T 3 , 5f u 6f 2 ,5f 3 ) T . 
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The dynamic model for the arm-length phenomenological variables is simply as follows 


d_ 
d t 
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(27) 


The dynamics for the clock variables are the same as Eq. [17 
The measurement matrix H\ 

■ "fe 

again explicitly write its first three rows as below 


_ is 18-by-14 and a bit different from that in the last section, we therefore 
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0...0 0 


/ carrier 
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v fi om ) 
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0 c c 0 


0 0 0 
0 0 0 


/ carrier r carrier 

3 _ J 1 
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/r rrier ^i 

/r m c 


o o 
0 0 


where we have omitted the step index k in the matrix components. 

We design a hybrid-extended Kalman filter with the phenomenological polynomial system model described above 
to process the simulated LISA measurement data. The results of the simulation are summarized in Fig. [8] Comparing 
these estimation errors with that given in the last section, we find that the overall performance of this phenomenological 
polynomial model is slightly better. Since this model is directly applicable to relativistic cases, we expect it to 
outperform the periodic model in general. 


histogram of estimate error of clock jitters histogram of arm length deviation histogram of frequency difference 



(a) 


(b) 


(c) 


FIG. 7: Histograms of the estimation errors of a hybrid-extended Kalman filter with a periodic system model in (a) 
differential clock errors (b) arm lengths and (c) differential frequency errors of the clocks for the laser link from S/C 

2 to S/C 1. 


histogram of estimate error of clock jitters 


histogram of arm length deviation 


histogram of frequency difference 



relative clock jitters [sec] 


arm length precision [m] 


frequency [Hz] 


(a) 


(b) 


(c) 


FIG. 8: Histograms of the estimation errors of a hybrid-extended Kalman filter with a phenomenological 
polynomial system model in (a) differential clock errors (b) arm lengths and (c) differential frequency errors of the 

clocks for the laser link from S/C 2 to S/C 1. 
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VIII. SUMMARY 

In this paper, we have designed and tested several alternative filtering dynamic models. The dynamic models with 
24-dimensional, 23-dimensional and 22-dimensional state vectors are comparable in reducing the measurement noise. 
The greatly simplified Kalman filter model with only clock variables performs a few times worse than the other filters 
as expected. All the filter models have successfully and significantly reduced the measurement noise. 

Most of the designed models use physical parameters, such as the positions and the velocities of the S/C, in the 
state vector to characterize the motion of the LISA constellation. The advantage of using these physical parameters 
is that they automatically fit into the physical laws, hence we can directly use Newton’s equations or relativistic 
equations as the dynamic equations to evolve the state vector of the Kalman filter. However, LISA interspacecraft 
measurements contain only relative quantities, such as the arm lengths, relative tangential velocities, relative clock 
jitter and relative frequency jitter. Therefore, the number of variables in the state vector is larger than the number of 
measurements, the physical variables in the state vector cannot be fully determined, and the measurement equations 
are very nonlinear. All of these add up complexities and numerical inaccuracies to the Kalman filter. 

Two phenomenological models have been designed and simulated. These models are superior to the models that 
directly use positions and velocities of the S/C, in the sense that they only use the observed relative quantities in their 
state vector. In turn, the state vector is much smaller, and the propagation matrices are much better conditioned. 
These models are also advantageous in reducing the nonlinearity of the dynamic equations and readiness to be extended 
to more complicated cases, such as full relativistic laser links. 


Appendix A: A limitation on the common clock drift 

As mentioned before, the differential clock errors can be estimated very accurately, whereas the common mode of 
the clock errors cannot be determined. This would result in errors in the time labels of the measurements. The errors 
in the time labels will introduce modulations to GW signals, hence they may affect the detection and the parameter 
estimation of GW signals. In this appendix, we try to estimate this effect and to set a limit on the permissible common 
clock drifts. 

We define a quantity cf t {t ) to characterize the timing stability 


a T (r) = ^((5T(k + l)-6T(k))*), (Al) 

where r is the nominal time increment between the sample k and the sample k + 1. Notice that this is different from 
the Allan variance [23] 


°i(j) 


1 / f S_fjk + 1) 

2 \ V / nom 


sm 


T((< 5T(k + 2) - 2 5T(k + 1) + 5T(k)) 2 ), 


(A2) 


which characterizes the frequency stability of the clock. 

For a GW signal with frequency /gw and a total observation time T 0 b s , if the mismatch caused by the common 
mode of the clock errors is less than e cycle, the effect on physical data analysis is negligible. Thus, we have a 
limitation on the common clock drift 


( j ~)/gw F 0 bs < e. (A3) 

For physical data analysis concern, a sampling time of r = 1 s will suffice. By considering the worst scenario, we take 
/gw = 0.1 Hz, T 0 bs = 10 8 s and e = 0.1, hence the timing stability should satisfy or(l) < 10 -8 . 
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